Anisotropic finite-size scaling analysis of a three-dimensional driven-diffusive system 
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We study the standard three-dimensional driven diffusive system on a simple cubic lattice where 
particle jumps along a given lattice direction are biased by an infinitely strong field, while those along 
other directions follow the usual Kawasaki dynamics. Our goal is to determine which of the several 
existing theories for critical behavior is valid. We analyze finite-size scaling properties using a range 
of system shapes and sizes far exceeding previous studies. Four different analytic predictions are 
tested against the numerical data. Binder and Wang's prediction does not fit the data well. Among 
the two slightly different versions of Leung, the one including the effects of a dangerous irrelevant 
variable appears to be better. Recently proposed isotropic finite-size scaling is inconsistent with our 
data from cubic systems, where systematic deviations are found, especially in scaling at the critical 
temperature. 

PACS number(s); 05.50.-|-q, 64.60.C, 05.70. Jk. 



I. INTRODUCTION 

Over the last decade or so, there have been many inves- 
tigations on nonequilibrium systems. One of the most of- 
ten studied is the driven diffusive system (DDS) for its 
simplicity of formulation and richness of novel properties. 
In particular, it is one of the few simple nonequilibrium 
systems showing phase transitions similar to that of an 
equilibrium statistical-mechanical system. The distinc- 
tion between DDS and an equilibrium system is a subtle 
one. Although both reach a stationary state of their re- 
spective stochastic dynamics in the long time limit, the 
DDS is generically nonequilibrium, defined only through 
its dynamics, while an equilibrium system can be alter- 
natively characterized by a Hamiltonian independent of 
the choice of dynamics. The standard DDS model is a 
lattice gas model governed by Kawasaki dynamics with 
a driving field. The driving field biases the motion of 
the particles in one preferred direction, so that under 
periodic boundary conditions it gives rise to a current 
along that direction. The existence of a steady current 
is a manifestation of the nonequilibrium nature of the 
stationary state. One physical realization of DDS is the 
superionic conductor. Certain flow properties of binary 
liquids under gravity, though more complicated due to 
hydrodynamic modes, are also similar to those of DDS. 

In this paper we focus our attention on the critical be- 
havior of the standard DDS. Mean-field theories |^ can 
give qualitative predictions of the phase transitions. But 
one of the achievements of the field-theoretic treatment of 
DDS is the exact determination of the set of critical 
exponents for all dimensions d > 2. After some contro- 
versies, the predictions in two dimensions have been con- 
firmed by extensive Monte Carlo simulations |^,||. The 
most difficult aspect of such tests is the fact that two 
length scales are involved — the correlation lengths (^ii 



and diverge differently in direction parallel to the 
field and in direction perpendicular to the field. Thus, 
one has to deal with anisotropic finite-size scaling with 
systems of various geometries. 

From a field-theoretic point of view, the d — 2 model 
is more complicated than in higher dimensions. This is 
because the usual 0^ coupling constant, denoted by u, is 
a dangerous irrelevant variable for d > 2, but it becomes 
marginal in ci = 2. While scaling arguments can predict 
the effect of u on finite-size scaling for d > 2 ||], there 
is little clue as to the existence and possible form of the 
associated logarithmic corrections in d = 2. Perhaps this 
explains why the agreements between previous tests 
and theory in d = 2 are not impeccable; small deviations 
from scaling may be due to the presence of small loga- 
rithmic corrections. For this reason, a more stringent yet 
practical test lies iw d — Z. 

Recently, not only the confirmations in d = 2 |5|j^ but 
also the validity of the field-theoretic approach itself have 
been questioned by Marro et al [Q. Isotropic finite-size 
scaling involving one scaling length was advocated. It is 
therefore our aim here to try to answer those questions 
by conducting a comprehensive test of the field-theoretic 
predictions in d = 3. Being free from the doubt of possi- 
ble logarithmic correction, consistencies between simula- 
tions and field theory would lend strong support to the 
latter. 

There are very old Monte Carlo simulation results 
of the three-dimensional DDS. One of the discoveries by 
Monte Carlo as well as by analytic works is the power- 
law long-range correlation of the particle density even 
in the high-temperature disordered phase Jl|,p|Jl^. This 
is shown as the manifestation of the violation of the 
fiuctuation-dissipation theorem in DDS. But the question 
of the validity of the field-theoretic results is far from 
answered. Monte Carlo simulation of the DDS model 
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is difficult for its long relaxation times and anisotropic 
correlations, due to local conservation and the external 
drive, respectively. In this article, we report a fairly ex- 
tensive Monte Carlo study using anisotropic finite-size 
scaling analysis similar to that in the two-dimensional 
case. A dominant feature of the anisotropics is the ap- 
pearance of an extra scaling variable, the "aspect ratio" 
S = l||^^/Lj_, with A = The exponents i/j^ and 

^"11 are associated with the correlation lengths and ^||. 
Simple data collapse among different samples is possible 
only when they have the same value of S, then the ordi- 
nary finite-size scaling forms are valid. There have been 
different predictions for the value of A: In d = 3, field 
theory predicts that A = 8/3, whereas Binder and Wang 
PI obtained A = 4. In this work, we mainly concentrate 
on samples with one fixed S, using the value 8/3 for A. 
The results support the field-theoretic predictions. Due 
to the huge demand on computer power, we could not 
study extensively the dependence on S, and the scaling 
at Tc, but only have some limited checks. We also test 
the assumption of isotropic scaling which corresponds to 
A = 1. Inconsistencies are found. 



II. THE DDS MODEL AND SIMULATION 
TECHNIQUE 

The model is defined on a simple cubic, fully periodic 
lattice of size L|| x L^. Each site on the lattice has a spin 
CTi = ±1. Equivalently, we can also consider the system 
as a lattice gas with local occupation variables. The total 
magnetization is set exactly at zero, and we assume ferro- 
magnetic interaction among nearest neighbors only, with 
a coupling constant J > 0. Equivalently, the particle oc- 
cupation is half filled and they attract each others. The 
system evolves according to the standard Kawasaki dy- 
namics of spin exchanges except with an extra ingredient 
due to an external "electric" field. We associate a posi- 
tively charged particle with an up spin and a hole with a 
down spin. Particle hoppings along the electric field are 
favored. Hereafter we will use the spin language. In sim- 
ulation, we consider only the extreme case of infinitely 
strong field, chosen to be in the +x direction. Thus, ex- 
change is always performed if we have a — ) pair along 
the +x direction and the exchange for (— ,+) is forbid- 
den. When the exchanges are perpendicular to the field 
(which happens 2/3 of the time), the field does not play 
any role. In that case, we compute the change in energy 
AE due to the exchange, accepting it with the Metropo- 
lis rate min (l,exp(-A£;/fcr)). From now on, we will set 
J/k to one. 

Since the model evolves very slowly via a conservative 
dynamics, a fast algorithm is of the utmost importance 
for the present undertaking. We used a multi-spin coded 
program to simulate 8 or 16 systems simultaneously, de- 
pending on the word length of the given machine. The 
method is similar to that of Kawashima et al , capa- 



ble of achieving a speed of about 0.1 /isec/(spin flip) with 
a typical workstation. 

We define the order parameter as 
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2L L^'^^IH' VmiM' + mO,lW, (1) 



where 



a{l, m, n) 



iy-iLi-iLi-i 

E E E' 

x=0 y=0 2=0 



e 



(2) 



The normalization is chosen such that (p — 1 for a slab 
geometry (the completely phase-separated configuration 
in the limit T 0). The following quantities are cal- 
culated: (i) the averaged order parameter m = ((/>), (ii) 
the "susceptibility" proportional to the fluctuation of the 
order parameter, 



LiiL 



X 



rsin(7r/Li) 

the susceptibility above the critical temperature, 



X 



Tsin(7r/i^) 



and (iii) the fourth-order cumulant. 
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(3) 



(4) 



(5) 



Note that g goes from 0.5 to I as temperature T goes 
from oo to 0. We will not report the results on x'] it 
yields no additional information as {(f)'^) appears to scale 
like (0)2. 

The computations are performed on a variety of work- 
stations, such as Alpha-stations, Pentium clusters, IBM 
SP2, etc. Our main results are obtained from a set of 
system sizes (Lj|,i_L) with (113,18), (193,22), (367,28), 

and (524,32), chosen such that S = L^^^^^/L^ is very 
close to a constant, ranging between 0.32703 to 0.32709. 
A second set with (59,35), (102,43), (122,46), (161,51), 
and (248, 60) are used mainly to confirm the result for 
Tc- The value of S for this set varies from 0.13171 to 
0.13182. Much more elongated geometries are also used 
to investigate the L[| oo behavior. A third set with 
cubic geometry, L = = L[| = 20, 30, 40, 50, 60, is 
used to test isotropic scaling. The lengths of runs are 10^ 
to 10* Monte Carlo steps per temperature, the longer for 
T closer to T^. We monitor the results until the system is 
well equilibrated before actually taking data. The total 
amount of CPU time spent is of the order of seven years 
on one IBM SP2 node. This gives an idea of the com- 
putational demand in achieving good statistics for this 
model. 
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III. ANISOTROPIC FINITE-SIZE SCALING 



xiT) 



AL 



(12) 



There have been two competing theories on the 
anisotropic finite-size scahng of the DDS. The first is that 
of Binder and Wang ||ll[]. Their results are based on a 
generalization of the one-dimensional Ginzburg-Landau 
type effective Hamiltonian for the very elongated geom- 
etry. The second is due to Leung based on the field- 
theoretic formulation 

Encouraged by the finite-size scaling of the standard 
Ising model above its upper critical dimension and the 
finite-size scaling at a Lifshitz point, Binder and Wang 
PI speculated on a scaling form for the driven diffusive 
model. The starting point is the assumption of an ef- 
fective functional for the local order parameter ^' in a 
quasi-one-dimensional geometry (L^ <C i||). 



dz 



1 /d"^' 

2 V dz>^ 
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where the exponent k"'^ = 2 + (5 — d)/3 which character- 
izes the singular term is introduced in such a way that 
the correlation-length exponent as obtained by field 
theory is reproduced at the Gaussian level. This leads to 
the scaling forms for the susceptibility and magnetization 
in three dimensions lOl at T = Tr 



m[T) = L-''^m{Ll'^/L^,tL\'^) 



g{T)^~g{Ll'^/L^,tL 



.3/4^ 



(13) 



(14) 



The above results imply that the thermodynamic limit 
has to be taken carefully. The field-theoretic results are 
understood to correspond to the case where L|| oc — > 
oo. Besides this limit, the quasi-one-dimensional limit 
L|| — > oo with L_|_ held finite is also of interest (we assume 
that there exists a unique value of Tc for these different 
ways of taking the thermodynamic limitj Since x does 
not depend on L|| in this case, from Eqs. (Q), (^) and (p^), 
we have three different predictions for the susceptibility 
at Tr 



Binder & Wang x~ Lj 
Leung (with u) L 
Leung (without u) x^ 



7/3 



(15) 
(16) 
(17) 



We will use Eqs. (0)-(^^ to check which version of the 
predictions for scaling describes the data of computer 
simulation better. 



x{Tc) 



(7) 



IV. ANISOTROPIC SCALING RESULTS 



m(Tc) 



L-'/'miLl/'/L, 



)• 



(8) 



On the other hand, Leung deduced the off- Tc finite- 
size scaling forms by generalizing the exact field-theoretic 
results for infinite system sizes to finite sizes The 
derivation was based on a combination of the renormal- 
ization group argument which treats 1/L± and l/i|| 
as two independent relevant variables, and scaling ar- 
guments which assume multiplicative singularities in the 



limit u L 



-2e/x 



0. 



Here 9 is proportional to the anoma- 
lous dimension of u. The fact that 6' = 1 — (5 — d)/3 > 
for d > 2 prescribes the dangerous nature of u above two 
dimensions. The main results read H 



xiT)^Lj/'x{L'/'/L^,tL'''), 



7/8x 



m(r) = L^'/'^m(ifVi±=^4^^), 



g{T)^~g{Ll^'/L^,tL]/'), 



(9) 



(10) 



(11) 



where t — [T — Tc)/Tc- If the contribution from the 
dangerous irrelevant variable u was ignored, we would 
have instead 



An accurate determination of the critical temperature 
Tc is important for a quantitative analysis of the critical 
behavior. We determine by the finite-size effect of 
the location Tpoak of the susceptibility peak in x- The 
data near peaks are fitted with a parabola to derive their 
heights and locations. First we consider the predictions 
in Eqs. (^ and (p^. Following well-known argument, 
the shift of Tr due to finite sizes is 



Tpcak(i|i, -S") = Tc + a(S')Tjj 



(18) 



where a{S) is a scaling function, h = 7/8 or 3/4. When 
S ~ T|^^/L^ is fixed, we have the usual shift of the 

peak locations. Fig. |^(a) shows Tpoak against 
for two sets of data with S = 0.327 (first set, lower 
part) and S = 0.1317 (second set, upper part), respec- 
tively. Least-square fit extrapolates to critical tempera- 
tures 4.859 ± 0.005 (first set) and 4.869 ± 0.005 (second 
set). If the exponent 3/4 is used in Eq. (|l8|), accord- 
ing to Eq. (|l|), the estimate shifts to higher values of 
4.870 ±0.005 (first set) 4.873 ±0.005 (second set). Thus, 
while the extrapolated Tc from the two sets of data agree 
within errors for both versions, the one with the exponent 
3/4 without the dangerous irrelevant variable correction 
appears to be marginally more consistent. The consis- 
tency in Tc for two different sets of data with fixed S 
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is a significant confirmation that A = 8/3 is the correct 
anisotropic scahng exponent. 

The peak heights are additional information which we 
can use. According to Eq. (||) and (|l^), the peak height 

scales with system size as Xmax ~ ^||^^ Xmax ~ L^\^^ 
if the variable u is or is not taken into account. The 
nice feature of the susceptibility maximum scaling is that 
it does not depend on the choice of the second scaling 
variable tL^^ in the scaling functions. In Fig. |2[ we plot 
the height vs. L|| in logarithmic scale. The insert shows 
shift of the susceptibility peaks. Very nice linear behavior 
with a slope of 0.86 ± 0.02 is obtained. This is a solid 
confirmation of Eq. (||) . 

For fixed S, the anisotropic finite-size scaling has the 
same form as the isotropic scaling, involving a prefactor 
in L and just one scaling variable of the form tL'^. We 
first look at the scaling of the fourth order cumulant. 
This quantity has a simple scaling since there is no pref- 
actor. Fig. 0(a) is the cumulant g plotted against scaling 

variable ^-^|| • If the dangerous irrelevant variable u was 
ignored, the exponent c would be 3/4. The correspond- 
ing scaling plot is Fig. |^(b). Both of them yield equally 
well data collapsing if different Tc are used. 

The finite-size scaling of the order parameter is pre- 
sented in Fig. IJ. Excellent scaling is found there when u 
is taken into account, using Eq. (|lO|). The upper branch 
is for t < and the lower branch is for t > 0. Fig. 3(a) 
shows that the upper branch has a slope of about 0.5 for 

7 /S 

large value of tL^^ , which is consistent with the expo- 
nent (3 — 1/2. If u is ignored, using Eq. ([l3| ) instead, 
Fig. 3(b) shows that the data collapse is not as good 
below Tc- 

The susceptibility data are shown in Fig. ^ in scaling 
form. The upper curve with a peak is for t < and 
the lower curve is for t > 0. The scaling is not as good 
as that for the order parameter at low temperatures. A 
plausible explanation for such deviations from scaling is 
that the low-T data may fall outside the critical region. 
If this is the case, the size of the critical region could then 
be estimated to be about 10% of Tc- This interpretation 
is supported by the same kind of deviation in the low-T 
tails in x ^-nd m in equilibrium Ising model jist , where 
the exponents are known exactly and thus cannot be the 
source of deviations. Similar behavior is also found in 
two-dimensional DDS [||. 

We now turn our attention to finite-size scaling at Tc, 
where only one scaling variable of the form L^^^^ /Lj_ is 
left. We simulate a wide range of system sizes and shapes, 
no longer restricting to ones with fixed Tj|/^/L_L. The ex- 
ponent A can be determined when ensemble averages for 
different (T||, Lj_) fall on one curve when plotted against 

t||^^/Lj_. Since large systems take very long time to 
equilibrate at Tc, we do not have very precise data. Nev- 
ertheless, it is sufficient to distinguish among alternative 
predicted scaling forms. First, the fourth-order cumulant 



at Tc = 4.860 is presented in Fig. ^. The prediction of Le- 
ung with or without u term has the same scaling variable 
i|^^/Tj^, while that of Binder and Wang is L^^^^^/L±. 
Isotropic scaling variable is T||/Lj^. The three sets of 
curves correspond to these three cases. Leung's scaling 
appears better. This plot also shows that for S w 0.14 
we get a maximum value in g. This value corresponds to 
geometries where two correlation lengths have the same 
ratio to the respective linear dimensions. For the as- 
sumption of T||/Lj_ as scaling variable, the data clearly 
do not scale well. This appears to be a strong evidence 
in favor of anisotropic scaling. The same type of plots 
assuming Tc — 4.872 does not give good scaling for all 
the choices of the scaling variables. This seems to im- 
ply that Tc — 4.860 is a better estimate for the critical 
temperature. 

In Fig. 1^ we show the scaling of magnetization at 
Tc = 4.860. Fig. 0(a) uses Binder and Wang scaling, 
(b) Leung's scaling with u correction, and (c) without u 
correction. Case (a) has large deviations; case (b) gener- 
ally scales better except at large i^^^/T_L. The scaling 
assuming Tc = 4.872 (not shown) looks worse for case 

(a) , and better for (c) at large scaling variable. Case 

(b) scales the best in both temperatures. This shows 
that the relative quality of scaling is not sensitive to the 
choice of Tc within its extrapolated range. The scaling 
function ifi{x) has the following asymptotic behaviors at 
large T|| or L± limit due to the sum of magnetization 
of totally independent regions. For small x, m oc l/i±, 
this implies m(x) '~ x for all three cases, for large x, we 
have m oc 1/T||^^, thus m{x) ^ x^ , y = —1/2, —1/6, 0, 
respectively for case (a), (b), and (c). Unfortunately, all 
three plots are more or less consistent with this asymp- 
totic slopes and thus alone it cannot give a sensitive test. 

Fig. ^show the corresponding scaling plots for the sus- 
ceptibility at the same choice of Tc. The trend is the same 
as in m. Although all three cases satisfy Eqs. (p^-(p^) 

for large T||^'^/ix, case (b) is most consistent with the 
notion of a scaling function x{x)- 

Finally, in Fig. we do a separate test of the pre- 
dictions of Eqs. (p^)-(p7|) for the very long geometry 
(i|l oo, L± finite). The critical temperature is taken 
to be Tc « 4.860. The T|| ^ oo limit is obtained by sys- 
tems (1280,20), (960,30), (320,40) together with smaller 
systems and extrapolated to large L||, assuming a 1/T|| 
convergence or similar power. The least-square fit to the 
data in Fig. ^ gives exponent 2.27 ± 0.07, in good agree- 
ment with Leung's scaling with u taken into account. 
However, if Tc — 4.872 is used, the exponent reduces 
to 2.04 ± 0.09. Thus the distinction with or without u- 
correction is not clear cut. On the other hand. Binder 
and Wang's scaling requires the exponent to be 3. This 
seems unlikely to be satisfied, due to the lack of data 
collapse in Fig. ^(a). This demonstrates that to pass a 
consistency check at Tc requires both a good data col- 
lapse and the correct asymptotic behaviors. 
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V. ISOTROPIC SCALING 

There have been arguments that DDS under 

infinitely large driving field could follow the normal 
isotropic finite-size scaling with one single correlation 
length exponent, although no specific prediction of scal- 
ing exponents are given. For completeness, we test it by 
analyzing cubic systems with L — 20, 30, 40, 50, and 60, 
assuming the usual finite-size scaling: 



X{T) = L-'/'^xitL'/''), 
m(T) = L-P/^mitL^/"), 



(19) 
(20) 
(21) 



A second scaling variable L||/L_l, equal to one in our 
data, is not explicitly written. In Fig. ^ we plot the lo- 
cation Tpcak of the susceptibility peak versus L~^/^ with 
two choices of v. Values of v outside this range result 
in unacceptable, nonlinear behavior. This plot therefore 
gives us an idea of the bounds of Tc. Peak extrapola- 
tion depends on the assumption on varying from 4.86 
to 4.90 for V — 0.67 to 1. The curve is not quite linear 
for any choice of v] this is expected if the true behavior 
is given by Eq. (|l^). The insert shows the intersections 
of the fourth-order cumulant. The intersections shift to- 
wards higher values from 4.83 to 4.86 as L increases. If 
we trust the larger system sizes, the estimate of Tc would 
be 4.86. This is the value barely in agreement with the 
extrapolation. 

Having determined T^^ the best off-Tc scaling is seen 
in m, achieved by choosing — 4.863, 1/v ^ 1.53, and 
P/v = 0.638 (see Fig. ^l]). However, this set of parame- 
ters is not the best choice for g. Figure 12 shows that with 
the same Tc and v as for to. Due to the fact that g curves 
are not intersecting at a unique T, systematic deviations 
are observed no matter what values of Tc and v are used. 
Next, we determine ^/v = 1.81 from the susceptibility 
peak height vs. system size. Relatively large deviation 
also occurs for the off-Tc susceptibility data (Fig. |l3|) . 

Even though the data more or less scale in tL^^" with 
a suitable set of exponents, we see systematic devia- 
tions and slight inconsistency among different quantities. 
The quality of data collapsing is inferior to that of the 
anisotropic scaling except perhaps for to. Since the crit- 
ical temperature estimated from isotropic scaling agrees 
with that of anisotropic scaling, the same data obtained 
for the purpose of testing anisotropic scaling at Tc can 
also be used to test isotropic scaling. The results do 
not support the latter. We have already seen that the 
fourth order cumulant g does not scale (Fig. ^, a ~ 1). 
In Fig. |lj we show mL^/" vs. T||/Tj^ at Tc. If the 
system is isotropic, it should scale. But clearly, we see 
large systematic deviations. Finally, in the quasi-one- 
dimensional limit, the equation analogous to Eqs. (p^)- 

( p^ ) is X c« Li^/^ = T^^, but the actual exponent from 
the data is much larger (Fig. ^). 



The reason for the above good fits off Tc and the poor 
fits at Tc is a consequence of our fitting the former first. 
This warrants the crossing at L^\/L± = 1 in Fig. |l^. The 
fact that there is no collapse elsewhere along the L^\/L± 
axis is a clear indication of the failure of the assumption 
of isotropic scaling. Thus, the moral of this exercise is 
that an apparently good fit from a partial test may be 
dangerously misleading in the case of scaling with two 
variables. 



VI. CONCLUSION 

We have performed a large-scale simulation for the 
standard driven diffusive model. We test the theoretical 
predictions of finite-size scaling against numerical data, 
as carefully, critically and completely as we can. The 
task is to decide which of the various competing theories 
is the most consistent with the data. There have been 
four theoretical proposals, characterized by different sets 
of finite-size scaling exponents. Two kinds of scaling tests 
are performed: at and off the critical temperature. Both 
are necessary because there are two scaling variables in- 
volved. The two versions due to Leung, with and without 
the dangerous irrelevant variable correction, fit the off-Tc 
data almost equally well if Tc is adjusted accordingly. As 
to the data taken at Tc, while they do not have enough 
quality for us to perform a stringent test, the scaling be- 
havior and susceptibility peak heights clearly favor the 
one including the dangerous irrelevant variable. Such 
tests also indicate that Binder and Wang scaling may 
not be valid. Regarding the recent proposal based on 
isotropic scaling, we find that the data can be cast in 
scaling forms with some effective exponents, but its the- 
oretical basis is more conjectural, and the quality of scal- 
ing is poorer especially at Tc. Thus, taking into account 
all aspects of scaling, we conclude that the field-theoretic 
prediction based on spatial anisotropies with dangerous 
irrelevant variable corrections is the most satisfactory. 
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FIG. 1. The susceptibility peak location rpcak(i||,S) as a 
function of L'^^; (a) b = 7/8, (b) b = 3/4. The system sizes 
(L||,Lx) are (113,18), (193,22), (367,28), (524,32) [solid cir- 
cle]; and (102,43), (122,46), (161,51), (248,60) [open square]. 
The limiting value Tpca,k{L^^ ^ (X,S) is an estimate of Tc- 




FIG. 2. Scaling of the susceptibility maxima. The slope is 
as predicted in Eq. (^, within error from the best fit 0.86. The 
insert shows the shift of the susceptibility peaks for {L\\,L±) 
being (113,18), (193,22), (367,28), and (524,32) from left to 
right. 
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FIG. 3. The finite-size scaling of the fourth order cumu- 
lant g with the scaling variable (a) tL^^ and Tc = 4.860; 

and (b) fLjJ/" and Tc = 4.872. The system sizes (L||,I/x) 
are (113,18) [diamond], (193,22) [triangle], (367,28) [square], 
(524,32) [circle]. 



FIG. 4. The finite-size scaling of the order parameter m 
with the scaling variable (a) tLy^ and Tc = 4.860; and (b) 

tL^^ and Tc = 4.872. The system sizes are the same as in 
the previous figure. 
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FIG. 5. The finite-size scaling of the susceptibility x with 
the scaling variable (a) tL^^ and Tc = 4.860; and (b) tLV'^ 
and Tc = 4.872. The system sizes are the same as in the 
previous figure. 
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FIG. 6. The scaling of the fourth-order cumulant at 
Tc = 4.860. The left sot of curves with a = 1/4 refers to the 
prediction of Binder and Wang; the middle set with a = 3/8 
is that of Leung; the right set assumes isotropic scaling a = 1. 
Same symbol means the same transverse system size L± =40, 
30 and 20, from left to right. 



9 




FIG. 7. The finite-size scaling of magnetization at 
Tc — 4.86 for (a) Binder and Wang scaling; (b) Leung scaling 
with u; (c) Leung scaling without u. Same symbol means the 
same transverse system size L± — 60, 40, 30, 20 and 10, from 
left to right. 



FIG. 8. The finite-size scaling of susceptibility at Tc = 4.86 
for (a) Binder and Wang scaling; (b) Leung scaling with u; 
(c) Leung scaling without u. Same symbol means the same 
transverse system size L± = 60, 40, 30, 20 and 10, from left 
to right. 
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FIG. 9. The susceptibility x a-t Tc in the limit — > cxa as 
a function of the perpendicular dimension Lx in logarithmic 
scales. The critical temperature used is Tc — 4.860 for the 
circles and Tc — 4.872 for open squares. 



FIG. 11. Scaling of the magnetization for the cubic sys- 
tems. The linear sizes are 20 (solid circle), 30 (square), 40 
(diamond), 50 (open triangle), 60 (solid triangle). 
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FIG. 13. Scaling of the susceptibility for the cubic systems. 
The system sizes are the same as in the previous figure. 




FIG. 14. Isotropic scaling for magnetization at Tc, using 
the same set of original data as in Fig. ^ 
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